Electromagnetic energy-momentum equation without tensors: 
a geometric algebra approach 



oo 

O 

o 



0^ 



Oh. 
I . 

o ■ 

Ph. 



> 

(N 

00 

en 

o 

oo 

O 



% 



Quirino M. Sugon Jr.* and Daniel J. McNamara 

Ateneo de Manila University, Department of Physics, Loyola Heights, Quezon City, Philippines 1108 
*Also at Manila Observatory, Upper Atmosphere Division, Ateneo de Manila University Campus 

e-mail: qsugon@observatory.ph 

July 10, 2008 



Abstract. In this paper, we define energy- momentum den- 
sity as a product of the complex vector electromagnetic field 
and its complex conjugate. We derive an equation for the 
spacetime derivative of the energy-momentum density. We 
show that the scalar and vector parts of this equation are 
the difi^erential conservation laws for energy and momentum, 
and the imaginary vector part is a relation for the curl of the 
Poynting vector. We can show that the spacetime derivative 
of this energy-momentum equation is a wave equation. Our 
formalism is Dirac-Pauli-Hestenes algebra in the framework 
of Clifi^ord (Geometric) algebra CUfi. 



1 Introduction 

The conservation laws for electromagnetic energy and mo- 
mentum are given in Simmons and Guttmann jj as 



E J = 



f- 



(1) 



1 as 



pE+jxB = -Vf/-— — --Heo(V-E)E+(E- V)E 



dt 

-fMo(V-H)H+(H-V)H). 



(2) 



where U is the energy density and S is the Poynting vector. 
Our aim is to unify these two equations. 

The standard way to unify Eqs. ((T} and ^ is through 
tensors as given in Jackson: [2] 



^^.e"" = 



c 



(3) 



where da is the spacetime derivative operator, O"'^ is the sym- 
metric stress tensor, Ffsx is the electromagnetic field tensor, 
and Jx is the four-current density. 

Yet, it is not obvious that the tensor equation in Eq. (|3]) 
is equivalent to the vector equations in Eqs. ([1} and ([2]). To 
prove this equivalence is straightforward but tedious. What 
we need is a formulation that would enable us to extract the 
two conservation laws in a single step. 

To answer this problem, we propose geometric algebra. [3] |4] 
This algebra combines the imaginary number i with the dot 
and cross products of vectors in a single associative vector 
product :i5j 

ab = a • b + i(a X b), (4) 



which looks like the Pauli identity in Quantum Mechanics [5] 
but without the matrices. 

One use of the Pauli identity is in the unification of the 
four Maxwell's equations in Electrodynamics: [71 [S] 



(ia + V)(E + iB) =p-i. 
c c 



(5) 



Separating the scalar, vector, imaginary vector, and imagi- 
nary scalar parts of Eq. ([5]) yields Gauss's, Ampere's, Fara- 
day's, and magnetic flux continuity laws. (For comparison, 
tensor calculus 9 and differential forms [I^ can only reduce 
Maxwell's equations into two.) 

So we ask: is there an energy-momentum counterpart to 
the Majcwell's equation? 

This question was answered by VoldQT] by left-multiplying 
the Maxwell's equation by the reverse (|) or conjugate of the 
electromagnetic field F = E -|- iB, and adding the resulting 
equation with its reverse: 



F'^idt + V)F + F'^idt + V)F = F^'S + S'l'F, 



(6) 



where S — p — j/c and the overdots determine the direc- 
tion of differentiation. The reversion operator, just like the 
Hermitian conjugation operator in Quantum Mechanics [12], 
changes the order of factors and operators-even the direction 
of differentiation. 

Yet, it is still not obvious that the scalar and vector parts 
of Eq. ^ are indeed Eqs. ([T} and ©. Showing this result 
is not trivial, because Eq. (|6]) must first be converted into a 
form involving the spacetime derivative of the electromagnetic 
energy-momentum densitv[T5] 



iFFt^(E + iB)(E-iB) = i(|Er 



IBI- 



-ExB. (7) 



To perform this differentiation properly, we shall use a theo- 
rem in Jancewicz|14) for the spatial derivative of a product of 
two vectors: 

V(ab) = (Va)b-a(Vb) + 2(a- V)b. (8) 

(In Jancewicz, b is replaced by b = ib; we only factored out 
the i.) 

We shall divide this paper into four sections. The first 
is Introduction. In the second section, we shall review 
geometric algebra and calculus within the framework of 
Hestenes's spacetime algebra in spacetime split form via Clif- 
ford (Dirac) algebra C/4 n. [15lll6j In the third section, we shall 



revisit Maxwell's equation and use it to derive the Energy- 
Momentum equation. We shall show that the scalar and vec- 
tor parts of the latter are the two conservations laws, while 
the imaginary vector part is a relation for the curl of the 
Poynting vector. The fourth section is Conclusions. 

2 Clifford (Geometric) Analysis 
2.1 Geometric Algebra 

The Clifford (Dirac) algebra CUfi is generated by four vectors 
eo, ei, 62, and 63 that satisfies the orthonormality relation 



GmGi/ "T Gi^Gn 



'iSuv, 



(9) 



for fiyiy — 0, 1,2,3. That is, e^ = 1 and e^e^ = — e^e^. We 
shall refer to eo as the unit temporal vector and to ei, 62, 
and 63 as the three unit spatial vectors. 

One important subalgebra of CZ4,o is the Pauli algebra C/3,0. 
If a and b be two vectors spanned by the three unit spatial 
vectors in Chfi, then by the orthonormality axiom in Eq. ([S}, 
we can show that a and b satisfy the Pauli identity: 



ab = a ■ b + i(a x b). 



(10) 



where the imaginary scalar i = 616263. In words, the Pauli 
identity states that the geometric product ab of two vectors is 
equal to the sum of their dot product a- b and their imaginary 
cross product i(a x b). 

In general, every element A in Chfi is expressible as a clifi'or 
or a linear combination of a scalar, a vector, an imaginary 
vector, and an imaginary scalar: 



A = A0 + A1+ iA2 + iAs 



(11) 



The spatial inverse or the automorphic grade involution [17] 
of A is defined in terms of the unit temporal vector eo = ° : 

i° = °AK (12) 

So by the orthonormality axiom in Eq. ((9)1 , we have 

A'' =Ao-Ai+iA2-iA3. (13) 

Notice that the spatial inversion operator changes the sign 
only of vectors and imaginary scalars. This differs from the 
reversion operator which changes sign only of imaginary vec- 
tors and imaginary scalars. [15] 

2.2 Geometric Calculus 
2.2.1 Time Differentiation 

The time derivative of the product of two vector functions a 

and b is 

9 , , . 9a, 9b ,^ ,, 

_(ab) = -b + a-. (14) 

The scalar and imaginary vector part of this equation are 

9 , , , (9a , 9b ,, ^, 

-(a.b) = ^-b + a.-, (15) 

9 , , , 9a , 9b ,^^, 

-(axb) = _xb + ax-, (16) 



2.2.2 Space Differentiation 

The spatial derivative operator V is defined as 

^9 9 9 9 ,,„. 

V= — =61-— +62^^— +63^^— , 17 

or 0x1 0x2 ax'i 

where we used the identity e^^^ = 6fe. Because V is a vector 
operator, we may use the Pauli identity in Eq. (|10|l to write 



Va = V-a + i(Vxa), 
aV = a-V + i(axV). 



(18) 
(19) 



Notice that Va is a function, while aV is an operator. 

For the spatial derivative of the geometric product of two 
vector functions a and b, we use overdot notation [19]: 



where 



V(ab) = (Va)b + Vab, 



• • 9b 9b 9b 

Vab = 61a- h 62a- h 63a- — . 

OTl 0X2 0X3 



(20) 



(21) 



If we employ the orthonormality axiom in Eq. @, Eq. (|2ip 
becomes 



Vab 



(ai6i — 0262 — 0363)61 



9b 

9x1 



+ (-aiei + 0262 — 0363)62 



9b 

dx3 
db 



-|-(-ai6i - 0262 + 0363)63^ — . (22) 

Adding and subtracting (a ■ V)b and rearranging the terms, 
we get 

Vab = -aVb + 2(a- V)b. (23) 

Hence, 

V(ab) = (Va)b-a(Vb) + 2(a- V)b, (24) 

which is Jancewicz's theorem in Eq. ([8}. 

To verify the correctness of the product rule in Eq. (I24|l . 
we separate its vector and imaginary scalar parts: 

V(a- b) - V X (a X b) 

= (V ■ a)b - (V X a) X b - a(V ■ b) 

+ax (V X b) + 2(a- V)b, (25) 

iV-(axb) = i(V X a) -b-ia- (V X b). (26) 

The second equation is a familiar identity, while the first can 
be verified through two other known identities: 120' 

V(a-b) = (a- V)b + (b- V)a + ax (V X b) 

+b X (V X a), (27) 

Vx(axb) = a(V-b)-b(V-a) + (b- V)a 

-(a-V)b. (28) 

2.2.3 Spacetime Differentiation 

An event r° is defined by its position r and time t: 

f° = (ct-f r)° = °(ct-r) = °f^ (29) 

where c is the speed of light. The square of the event is the 
Minkowski metric: 



which are the familiar identities in vector calculus. 



..o ~o -00 ~t -~t 2,2 I i2 



(30) 



(In Hestenes's spacetime algebra, the event a; = f°, so that 
3^70 ~ f°° = f — ct + r and x'^ = f°f° .) 

Corresponding to the event f° is the event derivative oper- 
ator 

where we used the identities Gq ^ = eo and (e^eo) ^ — eoe^. 
We can show that the square of event derivative operator is 
the d'Alembertian operator for wave equations. [21) 



3 Classical Electrodynamics 

3.1 Maxwell's Equation 

In a linear and isotropic medium characterized by the speed 
of light c = 1/y/Jii and the radiation resistance (^ = Wfije, 
the Maxwell's equation may be written as 



dE 
df° 



dE 
df 



Cf 



where 



E = E + iCH, 
J = (pc + j). 



(32) 



(33) 
(34) 



In words, the Maxwell's equation states that the event deriva- 
tive d/dr° of the electromagnetic field E is proportional to 
the event curent density j° .'22] 

Factoring out the unit temporal vector to the right of Eq. 
((321) . we get 



That is. 



ld_ 
cdt 



dE .t 



-V)(E + iCH)=C(pc-j) 



(35) 



(36) 



Separating the scalar, vector, imaginary vector, and imagi- 
nary scalar parts of Eq. (|36|) . we obtain Maxwell's equations: 



(37) 
(38) 



V-E = Cpc, 

if -cv.H . a 

i(-^+VxE) = 0, 
iC(V • H) = 0. 

3.2 Energy-Momentum Equation 

Let us define the event momentum density as 



(39) 
(40) 



c^ c c 



-±eE&° 
2c 



- — eE°E. 
2c 



where 



U 



i(6|Er + M|H|^) 



S = Ex H 



(41) 

(42) 
(43) 



are the energy density and the Poynting vector, respectively. 
The square of Eq. (|4T|) is 






(44) 



If this quantity is zero, then U — |S|/c, or equivalently, 
£ = |p|c, which is the energy-momentum relation for light-like 
particles. 

The event derivative d/df° of the energy-momentum den- 
sity S°/c^ is 



1 dS° 
(? df° 



._Le°^(£#)°, 
2c df^ ' 



(45) 



This leaves us with the problem of differentiating EE'' . 

By the chain rules in Eqs. (|14|> and (|24|> . we can show that 
the temporal and spatial derivatives of EE^ are 



d_ 
dt 



{E&) 



dE-^, dE^ 



(46) 



V(i5£t) = {\/E)E'' -eve'' + 2{E-\/)e\ (47) 

(48) 



where 

£ • V = E • V + iCH ■ V 
Adding Eqs. (|46l) and (|47l) yields 

^(E&) = ^&+E^ + 2{E-\/)El (49) 
or or or' 

If we substitute Eq. (|49)l back to Eq. (|45} and use the 
Maxwell's equation in Eq. (|32|l . we get 

dS° 
df° 



l(rE+''Er)^-^{E-V)E, (50) 



after dividing by —2^ and using the spatial inversion property 
of the temporal vector. Juxtaposing the two unit temporal 
vectors to cancel each other, and taking the spatial inverse of 
the resulting equation, we arrive at 



as; 

dr 



-hj^& + Ej)-l{E-V)E\ 



(51) 



We shall refer to Eqs. (|50|l and (|51[) as the event-momentum 
and energy-momentum equations, respectively. 
Expanding Eq. ([51) , 



(i|+V)(ct/ + S) 



(pc-j)(-E + jCH) 
(E + iCH)(pc+j) 



--E- V(-E + iCH) 
-iH- V(-E-hiCH). (52) 

and separating the scalar, vector, imaginary vector, and imag- 
inary scalar parts, we get 

dU 



dt 



-l-V-S 



cVC/ + 



-J • E, (53) 

+ i(E-V)E + C(H-V)H 

-CJ X H, (54) 

J(VxS) = i(-^H-(E- V)H-h(H- V)E),(55) 



19S 
c 9^ 







-C(-j-H + H.j) 



(56) 



The first equation is the conservation of energy. The second 
is the conservation of momentum, which may be converted 
to that in Eq. ^ by adding Gauss's law and magnetic fiux 
continuity law. The third is a relation for the curl of the 
Poynting vector S, which follows from the product rule in 
Eq. ([28p. The fourth is insignificant. 



4 Conclusions 

In this paper, we used Hestenes spacetime algebra in space- 
time split form to write down Maxwell's equation, which 
states that the event derivative oaf the electromagnetic field is 
proportional to the event current density. We showed that the 
four Maxwell's equations may be extracted from this equation 
by employing the spatial inversion of the temporal vector and 
Pauli expansion for the geometric product of vectors. 

We defined the event momentum density in terms of the 
product of the electromagnetic field with its spatial inverse. 
We used Jancewicz's theorem for the differentiation of a prod- 
uct of two vectors together with the Maxwell's equation to de- 
rive the event momentum equation, whose scalar and vector 
parts are the conservation laws for energy and momentum, 
and whose imaginary vector part is a relation for the curl of 
the Poynting vector. We can show that the event derivative of 
the event momentum equation would yield its corresponding 
wave equation. 
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